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Abstract We outline the main features of nuclei acceleration at supernova remnant 
forward shocks, stressing the crucial role played by self-amplified magnetic fields 
in determining the energy spectrum observed in this class of sources. In particular, 
we show how the standard predictions of the non-linear theory of diffusive shock 
acceleration has to be completed with an additional ingredient, which we propose 
to be the enhanced velocity of the magnetic irregularities particles scatter against, 
to reconcile the theory of efficient particle acceleration with recent observations of 
gamma-ray bright supernova remnants. 



1 Introduction 

Supernova remnants (SNRs) have been regarded for many decades now as the 
sources of Galactic cosmic rays (CRs), both because of their energetics and of the 
fact that strong shocks are expected to naturally accelerate particles with power-law 
energy distributions, according to the so-called Fermi mechanism. 

Clear-cut evidences of particle acceleration in SNRs have been found many years 
ago in radio to X-ray synchrotron emission, therefore attesting the presence of elec- 
trons with energies as high as 10-100 TeV. On the other hand, direct evidences of 
proton acceleration are much less clear: the most prominent signature of hadronic 
acceleration has been individuated in the y-ray emission from the decay of neutral 
pions produced in nuclear interaction between relativistic nuclei and the interstellar 
medium fPHl . The first SNRs detected in the TeV y-rays, however, did not unravel 
the question, since in this energy range it is hard to disentangle whether the ob- 
served spectrum (typically exhibiting a cut-off around ~10 TeV) is the result either 
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of the hadronic mechanism described above or of inverse-Compton scattering (ICS) 
of relativistic electrons on some photon background. 

Coupling the TeV observations of Cherenkov telescopes as HESS, VERITAS 
and MAGIC with the GeV observations of the Fermi satellite, however, it has been 
possible to assess the nature of the emission, at least in some paradigmatic cases: 
SNR RX J 17 13.7-3946 turned out to be mostly consistent with a leptonic scenario 
H] while, on the contrary, Tycho's SNR has been convincingly modeled according 
to the hadronic hypothesis fl2l . The main element allowing such a distinction is the 
slope of the observed y-ray spectrum: for a oc E~ q CR energy distribution, ICS is 
expected to produce a flatter power-law photon spectrum, °c £ (? +1 )/ 2 ^ w hile pion 
decay and bremsstrahlung produce a spectrum parallel to the parents' one, « Ey q . 

Both GeV and TeV observations are — in basically all the cases — consistent 
with spectra steeper than E~ 2 , as showed in Fig. [TJ(see ref. ||4) for more details). 
The only SNRs showing a GeV slope smaller than 2 are in fact the aforementioned 
RX J1713. 7-3946 and Vela Jr., whose emission may be accounted for by invoking 
a strong photon background making ICS dominate over pion decay. When the spec- 
trum is steeper than E~ 2 across many energy decades, however, a hadronic mecha- 
nism is favored, instead, since bremsstrahlung is typically negligible in SNRs. 
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Fig. 1 Slopes q-f of the y-ray spectra as inferred in the GeV/TeV band for both confirmed and 
candidate SNRs, as in the legend. Data from table 1 in ref. (3). 
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But this behavior has another important consequence in terms of the process 
responsible for particle acceleration (of both electrons and protons, which share 
the same spectrum in rigidity): it has to be able to produce particle spectra steeper 
than E~ 2 . This apparently harmless requirement is at odds with any theory of effi- 
cient particle acceleration at SNRs developed in the last 30 years since, according 
to Fermi's mechanism, the expected slope is a function of the shock compression 
ratio only, and namely q = . For strong shocks r = = 4, with y = 5/3 the 
gas adiabatic index and, in turn, q = 2. When acceleration is efficient, though, one 
could expect the shock to be modified by the back-reaction of the accelerated parti- 
cles, whose pressure may become comparable with the ram one. The simplest way 
to account for this effect is considering that CRs are relativistic particles forming a 
gas with adiabatic index y cr = 4/3, which would provide a compression ratio r = 7 
and in turn q = 1.5. The escape of high-energy particles may in addition make the 
shock behave as partially radiative, increasing r even further and pushing q down 
to 1 — 1.2. In any case, a very solid prediction of the non-linear theory of diffusive 
shock acceleration (NLDSA) is that the CR spectrum has to be invariably flatter 
than E~ 2 at the highest energies, and being flatter and flatter for larger and larger 
acceleration efficiencies (for more details on CR modified shocks see, e.g., ref. iFFD ). 

Therefore, independently of the possible origin (leptonic or hadronic) of the ob- 
served GeV emission, Fermi-LAT's observations force us to rethink our theory of 
CR acceleration, challenging the paradigm which requires SNR shocks to channel 
fraction as large as 10-30% of their kinetic energy into accelerated particles in order 
to be the potential sources of Galactic CRs. 



2 Magnetic field amplification and particle acceleration 

A fundamental tile in the mosaic of the comprehension of particle acceleration in 
(young) SNRs has also come from the detection of narrow X-ray bright rims imme- 
diately downstream of their forward shocks. The non-thermal nature of their spectra 
has been accounted for as due to synchrotron emission of relativistic electrons, while 
their narrowness (typically less than O.Olpc) points to magnetic fields as large as a 
few hundreds jiG, almost two orders of magnitude larger than the typical interstellar 
one. Such an evidence has been welcomed by theorists for several reasons. 

• The super-Alfvenic streaming of CRs upstream of the shock has been predicted 
to lead to the excitation of several magnetic modes via plasma instabilities |[T4l l2l. 

• A higher level of magnetization enhances particle diffusion and allows the 
achievement of larger energies, arguably up to a few PeV, namely enough to 
account for the knee observed in the Galactic CRs detected at Earth 0. 

• The pressure in magnetic turbulence may be so large to overcome the gas one, 
therefore preventing the CR pressure to modify the shock too severely (8). 

• Finally, magnetic field self-generation may have a key role also in determining 
the CR spectral slope, and therefore in reconciling NLDSA theory with y-ray 
observations, which is what we want to outline here (also see ref. 0). 
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Without delving into the details of a kinetic theory of plasmas, we can state that, 
being the accelerated particles and the magnetic fields coupled through resonant 
scattering, the whole system has to gradually reduce the relative velocity of CRs, 
whose diffusion velocity is basically the shock velocity in the upstream reference 
frame, u, and of magnetic irregularities, whose phase velocity is of the order of the 
Alfven speed, va = Bj \J\np. In other words, magnetic field amplification might 
be considered as an inevitable consequence of an efficient CR acceleration while, 
on the other hand, an enhanced level of magnetization is a necessary condition for 
an efficient particle diffusion and acceleration: if the initial Alfvenic Mach number 
Ma = u /va is large (as it typically is is in the interstellar magnetic field), it tends to 
be reduced by a non-linear interplay between particles and fields. 

We model such an interplay by adopting the semi-analytic formalism for NDLS A 
outlined in ref. J3] and references therein, which proved to be quick but very ac- 
curate for non-relativistic shocks @. We solve the equations of hydrodynamics 
coupled with a kinetic description of the non-thermal particles (via the diffusion- 
convection equation, eq. 2.5 in ref. f5|) and a fluid description for the generation 
of magnetic turbulence via resonant streaming instability (eq. 2.18 in ref. Q). The 
dynamical feedback of both accelerated particles and amplified magnetic fields are 
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Fig. 2 Top panel: pressure in CRs at the shock, P cr , expressed as a fraction of the ram pressure 
pu 2 , and predicted downstream magnetic field, Bj, as a function of the injection efficiency 77. The 
environmental parameters are those relevant for a typical S NR in the early Sedov stage (see ref. (5] 
for more details). Bottom panel: CR spectral slope obtained when the Alfven velocity is calculated 
in the self-amplified magnetic field SB (solid line with triangles) or according to a standard NLDSA 
theory assuming M^(Bq) ;§> 1 (dot-dashed line). 
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self-consistently retained in this approach, as well as the escape of the highest en- 
ergy particles from the upstream boundary [0. 

The key role of magnetic field amplification in determining the CR spectral 
slope is to enhance the velocity of the scattering centers, i.e., the phase veloc- 
ity of the magnetic waves, ~ v A oc §B. Since the compression ratio shaping 
the CR spectra is not the one of the fluid, but rather the one of the veloc- 
ity of their scattering centers ahead and behind the shock, when becomes 
a non-negligible fraction of u, it is in principle possible for diffusive shock 
acceleration to produce spectra steeper than E~ 2 . 

More precisely, the compression ratio felt by the fluid between upstream and 
downstream is the usual r = u u /u e i > 4 for strong, CR modified shocks, while the 
compression ratio felt by CRs is instead r = (u u + va,u) / (ud + vaa)- According to 
the argument above, however, the coupling between CRs and magnetic field is ex- 
pected to produce waves streaming against the fluid in the upstream (in the direction 
opposite to the CR gradient), while downstream the turbulence should be efficiently 
isotropised {v^.d ~ 0). 

The net effect is that f = (1 — M^)r and, in turn, if Ma(8B) becomes small 
enough, the predicted slope of the CR spectrum q = 4±y may become larger than 
2, even for strong shocks and large CR pressures. On the other hand, a steeper 
CR spectrum allocates less pressure (which in turn generates less magnetic field), 
therefore we need a non-linear calculation to understand which effect dominates the 
determination of the spectral slope when acceleration is efficient: the pressure in 
CRs, which makes the shock more compressible, or rather the reduced jump in r. 

Since CR injection has not been understood from first principles, yet, we can- 
not but parametrize it: we allow 77, i.e., the fraction of particle crossing the shock 
with a momentum large enough to be injected into the acceleration mechanism, to 
vary between ~ 10~ 7 and ~ 0.1. Such a large range of injection efficiency reflects 
different regimes of CR acceleration, spanning from unmodified shocks (where the 
pressure in CRs, P cr is expected to be negligible with respect to the ram pressure), 
to strongly modified ones, in which the shock structure is deeply affected by CRs 
and in which magnetic field amplification is very effective. At saturation, in fact, we 
have that P B = 8B 2 /8n ~ P cr /M A (see §2.1 in ref. J3) for the exact dependence). 

Our results are showed in Fig. |2]for a SNR about 3000 yr old, i.e., in the early 
Sedov stage. For a much wider discussion of the SNR evolution in terms of its non- 
thermal content the reader may refer to refs. [0 |5]. In the top panel, P cr and the 
downstream field B e / are showed as function of 77 : they reproduce the test-particle 
case for 77 < 10~ 5 , show a linear increase for about a decade and finally saturate to 
about 10-20% and 100-200^iG, respectively, above 77 ~ 10~ 4 . 

The very reason of this saturation, which among other things assesses the mild 
dependence of our findings on 77, is that the whole system reacts to high injection ef- 
ficiencies by steepening the CR spectrum in order to accommodate almost the same 
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energy in non-thermal particles, as demonstrated in the bottom panel of Fig.|2](solid 
line). Spectra as steep as E~ 2 - 2 — £~ 2 5 correspond to the highest acceleration effi- 
ciencies and to the fields inferred in the downstream of young shell SNRs (see figure 
3 in ref. 0). The spectral slope obtained is also compared with the one predicted 
by a two-fluid theory of NLDSA, q = (3y e ff — 1 )/2, where the CR contribution is 
accounted for via a compression ratio modulated by an effective adiabatic index of 
the gas+CR fluid ifTOl , namely Jeff = fzfp 1 (dot-dashed line in the bottom panel 
of Fig. [2J. It is worth recalling that the latter slope, already clearly inconsistent with 
the observations reported in Fig. [TJ represents an upper limit to the one which would 
be predicted by taking into account also the escape of CRs, which would make the 
shock partially radiative, in principle increasing r well above 7 iFFD . 

The only way to accommodate the typical efficiencies (10-20%) required for 
SNRs to be the sources of Galactic CRs, and the high levels of magnetization 
observed (a few 100/iG), with CR spectra steeper than E~~ is to modify the 
standard theory of NLDSA. 

We propose here that accounting for the features of the magnetic field pro- 
duced by the CR streaming is a possible way to build a consistent picture of 
efficient NLDSA in SNRs, actually reversing the usual trend which relates 
larger and larger CR acceleration efficiencies with flatter and flatter spectra. 
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